arXiv:1508.03434vl [hep-ph] 14 Aug 2015 



ELSEVIER 


Available online at www.sciencedirect.com 


ScienceDirect 


Nuclear Physics B Proceedings Supplement 00 (2015) 1-^ 


Nuclear Physics B 
Proceedings 
Supplement 


Solving the Balitsky-Kovchegov equation at next to leading order accuracy 


T. Lappi‘‘’'’, H. Mantysaari^* 


^Department of Physics, University of Jyvdsky Id, P.O. Box 35, 40014 University of Jyvdskyld, Finland 
^Helsinki Institute of Physics, PO. Box 64, 00014 University of Helsinki, Finland 


Abstract 

We solve the Balitsky-Kovchegov small-.r evolution equation in coordinate space. We hnd that the solution to the 
equation is unstable when using an initial condition relevant for phenomenological applications at leading order. The 
problematic behavior is shown to be due to a large double logarithmic contribution. The same problem is found when 
the evolution of the “conformal dipole” is solved, even though the double logarithmic term is then absent from the 
evolution equation. 
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1. Introduction 


2. The BK equation at NLO 


The Color Glass Condensate [I] effective theory of 
QCD at high energy has been shown to be in good 
agreement with large amount of experimental data on, 
for example, deep inelastic scattering, single and dou¬ 
ble inclusive particle production and exclusive vector 
meson production, see e.g. 13 sa. There are two 
main ingredients in these calculations; first, one needs 
the small-x evolution equation such as the Balitsky- 
Kovchegov (BK) equation which describes the 

evolution of the dipole-target scattering amplitude as 
a function of energy, or equivalently, Bjorken-x. The 
dipole amplitude at initial Bjorken-x is the second nec¬ 
essary input, and it can not be obtained from perturba¬ 
tive calculations but it must be fit to experimental data. 
The fact that one can indeed obtain a good description 
of the precise combined HERA single inclusive data JT) 
has been one of the tightest tests for the CGC 13|8l. 

An important next step for the description of the sat¬ 
uration phenomena from the CGC framework is to de¬ 
velop the CGC theory to the next to leading order ac¬ 
curacy. First steps in this direction have been taken by 
deriving e.g. the photon impact factor a and single in¬ 
clusive cross section ifTOll at this order. The NLO BK 
equation is also known im , but no solution to it existed 
before our work II3. 


The BK evolution equation for the dipole operator 
S, which is a correlator of Wilson lines U such that 
S(x - y) - l/Nc{ Tr t/^(x)t/(y)). At NLO accuracy the 
equation reads 

dyS(r) = ® [5(A)5(T) - 5(r)] 

+ ^^^2 ® [S{X)S{z - z')Sir) -SiX)SiY)] 

+ ^^Kf®S(Y)[S(X')-S(X)]. (1) 

The quark and the antiquark of the parent dipole are at 
transverse positions x and y, and the daughter dipole 
sizes are A = |x - z|, T = \y - z\. A' = |x - z'\ and 
Y' - \y - z'\, and r is the size of the parent dipole. The 
convolutions ® are taken by integrating over the daugh¬ 
ter dipole sizes (z in Ki and both z and z' in K 2 and Kf). 

The kernel Ki includes the leading order BK kernel, 
the mnning coupling part and an aj correction, as 
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We implement the running coupling by replacing the 
terms proportional to the p function coefficient by the 
Balitsky running coupling prescription from Ref. d. 
The kernel Ki then reads 


asNc _ as(r)Nc 
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We will later refer to the part proportional to 
InXVr^ln Y^/r^ as the double logarithmic term. 

The kernels K 2 and Kf are combinations of rational 
expressions of transverse separations and a logarithm 
\nX^Y'^ liX'^Y^). Note that this logarithm vanishes in 
the small parent dipole limit where » y, in contrast to 
the double logarithm. The coupling constant cTs is eval¬ 
uated at the scale set by the parent dipole, as it is the 
only external scale. For explicit expressions, we refer 
the reader to Refs. iniiii. 

As an initial condition for the NLO BK equation we 
use a modified McLerran-Venugopalan (MV) model 


N{r) - 1 - 5 (r) = 1 - exp 
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(4) 

Here the MV model is modified by introducing an 
anomalous dimension y which controls the power-like 
tail of the dipole amplitude at small dipole sizes. The 
leading order fits to the HERA data prefer JS) values 
of y ~ 1 . 1 , which then reduces during the evolution to 
y ~ 0.8. The constant Qs.o parametrizes the saturation 
scale at initial Bjorken-x. In this work, we do not seek 
for parameter values that are compatible with the ex¬ 
perimental data. In practice, controls the relative 
importance of the NLO terms as it scales the value of 


3. Solution to the NLO BK 

In Fig. [TJwe show the evolution speed dyN{r)fN(r) 
as a function of the dipole size for the MV model 
( 7 = 1 ) initial condition. At small initial saturation 
scales 2 s,o/Aqcd, when the strong coupling constant 
and the NLO corrections are largest, the evolution speed 
is negative at all dipole sizes. With smaller values of 
(larger saturation scale) the evolution speed turns nega¬ 
tive at small dipole sizes when r « i/Qs- 



Figure 1: Evolution speed of the dipole amplitude at initial condition 
(MV model with 7 = 1) with different values for the initial saturation 
scale. 


The negative evolution speed is unphysical, as it cor¬ 
responds to having an unintegrated gluon distribution 
that decreases when it is probed at smaller and smaller 
X. However, having dyNIN ~ In r in the small r limit 
is a signal of mathematical instability, as in that case 
there is a small (but finite) r below which the dipole am¬ 
plitude becomes negative in one step dy of the rapidity 
evolution. On the other hand, the definition of the dipole 
amplitude N(x -y) - 1 - l/Nd Tr UHx)U(y)) requires 
that N(r) —> 0 in the limit r —> 0. Also, if the dipole 
amplitude does not satisfy this requirement the z inte¬ 
gral in the leading order equation does not converge. In 
our numerical analysis we impose by hand a constraint 
N(r) > 0. 

Let us then trace back the origin of the negative evo¬ 
lution speed. In Fig. we show contributions to dyNfN 
originating from the different terms of the NLO BK 
equation. We observe that the double logarithmic term, 
which is part of the kernel Ki, is the one that drives the 
evolution speed negative. The other NLO corrections 
also decrease the evolution speed but do not cause the 
problematic dyN jN ~ In r behavior. 

It has been argued in Ref. El that the double loga¬ 
rithmic contributions should be resummed to all orders. 
The resummation effectively removes the double loga¬ 
rithmic term from the kernel and multiplies the lead¬ 
ing order BK kernel P' jX^Y^ by an oscillatory factor, 
which expanded to order gives the double logarith¬ 
mic term to the kernel Ki. The initial condition is also 
modified by the resummation procedure. We implement 
this resummation in our analysis and show in Fig.[^the 
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Figure 2: Contributions to the evolution speed of the dipole amplitude 
originating from dilferent terms of the NLO BK equation at the initial 
condition (MV model with y = 1). 


Figure 3: Evolution speed of the dipole amplitude obtained by re¬ 
summing the double logarithmic terms to all orders. Result is shown 
at initial condition (MV model with y = 1), and the contributions from 
the resummation and from the other NLO terms are shown separately. 


evolution speed obtained by solving the NLO BK equa¬ 
tion with resummation. We find that when the double 
logarithmic corrections are resummed to all orders the 
evolution speed turns negative at significantly smaller 
dipoles, but the negativity problem is not completely 
solved. However, it was later shown in Ref. ifTSi that 
there is also a single log term that should be resummed. 

Let us then study the dependence on the anoma¬ 
lous dimension y. We define the dipole size dependent 
anomalous dimension as yir) - dlnN(r)/dlnr2 and 
show it in Fig. at different rapidities using different 
values for the parameter y in the initial condition. With 
y = 1 in the initial condition, the anomalous dimension 
increases rapidly in the evolution, which is a signal of 
the dipole amplitude approaching a step function form 
N(r) ~ 0(r - 1 / Qs) (recall that we force N(r) > 0). This 
unstable behavior is not seen with y = 0.6 in the initial 
condition within the rapidity range studied here. Indeed, 
with y = 0.6 the evolution speed actually remains pos¬ 
itive at small dipoles (rAqcD ^ IQ-® included in the 
numerical calculation) within the studied rapidity inter¬ 
val. 


This non-conformal term also appears in the evolution 
equation derived in A = 4 Super Yang-Mills theory, as 
shown in Ref. Ea. It is therefore interpreted as an arte¬ 
fact of a cutoff that breaks the conformal invariance. 

In Ref. OH the authors present a way to restore the 
conformal invariance by writing the evolution equation 
for a conformal dipole defined as 


= 5(r) - 


A-jC- 
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d^z 


In 
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x[S(X)S(Y)-S(r)]. 


( 5 ) 


Here a is a dimensionful constant that drops out of the 
evolution equation. Writing the evolution equation for 
the conformal dipole removes the double logarithmic 
term from Ki and introduces a new term proportional 
to Inr^ in K 2 . We have solved the NLO BK equation 
for the conformal dipole, and as a result find that quali¬ 
tatively the evolution does not differ from the evolution 
of the “non-conformal” dipole. As shown in Ref. 02l. 
the new In r^ logarithm now drives the evolution speed 
negative at small dipoles. 


4. Conformal dipole 

The Wilson lines are by dehnition conformally invari¬ 
ant, and thus one expects the evolution equation for the 
dipole operator S to also respect this symmetry if the 
running of the coupling is neglected. However, in the 
NLO BK equation the double logarithmic term, which is 
responsible for the unstable behavior, is not conformal. 


5. Discussion 

We have obtained the first numerical solution to the 
Balitsky-Kovchegov equation at next to leading order, 
and the NLO corrections are shown in decrease the evo¬ 
lution speed. This is expected, as the leading order fits 
to the deep inelastic scattering data seem to favor rela¬ 
tively small evolution speeds. The double logarithmic 
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Figure 4: Anomalous dimension y{r) at rapidities y = 1,5,30 with different values for parameter y in the initial condition. 


term in the NLO BK equation, however, causes an un¬ 
stable behavior at small dipoles when an initial condi¬ 
tion relevant to phenomenological applications (at least 
at leading order) is used. The resummation procedures 
shown in Refs. oiiisi have potential to improve the 
situation and bring the NLO BK equation closer to the 
point where it can be used in phenomenological appli¬ 
cations. In addition to the resummation schemes, also 
the role of the kinematical constraint to the BK equa¬ 
tion iniiisi should be better understood. 

Acknowledgements 

This work has been supported by the Academy of 
Finland, projects 267321 and 273464, the Graduate 
School of Particle and Nuclear Physics (H.M.) and by 
computing resources from CSC - IT Center for Science 
in Espoo, Finland. 

References 

[1] F. Gelis, E. lancu, J. Jalilian-Marian, R. Venugopalan, The Color 
Glass Condensate, Ann. Rev. Nucl. Pail;. Sci. 60 (2010) 463- 
489. arXiv:1002.0333 doi:10.1146/annurev.nucl. 
010909.083629 

[2] T. Lapp!, H. Mantysaari, Forward dihadron correlations in 
deuteron-gold collisions with the Gaussian approximation of 
JIMWLK, Nucl.Phys. A908 (2013) 51-72. arXiv: 1209.2853 
doi:10.1016/j.nuclphysa.2013.03.017 

[3] T. Lappi, H. Mantysaari, Single inclusive particle production at 
high energy from HERA data to proton-nucleus collisions, Phys. 
Rev. D88 (2013) 114020. arXiv: 1309.6963 doi: 10.1103/ 
PhysRevD.88.114020 

[4] T. Lappi, H. Mantysaari, J/'¥ production in ultraperipheral 
Pb+Pb and p+Pb collisions at LHC energies, Phys. Rev. 
C87 (2013) 032201. arXiv:1301.4095 doi:10.1103/ 
PhysRevC.87.032201 

[5] 1. Balitsky, Operator expansion for high-energy scattering, Nucl. 
Phys. B463 (1996) 99-160. arXiv:hep-ph/9509348 doi: 
10.1016/0550-3213(95)00638-9 


[6] Y. V. Kovchegov, Small-x F 2 structure function of a nu¬ 
cleus including multiple pomeron exchanges, Phys. Rev. D60 
(1999) 034008. arXiv:hep-ph/9901281 doi:10.1103/ 
PhysRevD. 60.034008 

[7] F. Aaron, et al.. Combined Measurement and QCD Analy¬ 
sis of the Inclusive e^p Scattering Cross Sections at HERA, 
JHEP 1001 (2010) 109. arXiv:0911.0884 doi:10.1007/ 
JHEPOl(2010)109 

[8] J. L. Albacete, N. Armesto, J. G. Milhano, P. Quiroga-Arias, 
C. A. Salgado, AAMQS: A non-linear QCD analysis of new 
HERA data at small-x including heavy quarks, Eur. Phys. J. 
C71 (2011) 1705. arXiv:1012.4408 doi:10.1140/epjc/ 
S10052-011-1705-3 

[9] G. Beuf, NLO corrections for the dipole factorization of DIS 
structure functions at low x, Phys. Rev. D85 (2012) 034039. 
arXiv: 1112.4501 doi : 10.1103/PhysRevD. 85.034039 

[10] G. A. Chirilli, B.-W. Xiao, F. Yuan, Inclusive Hadron Produc¬ 
tions in pA Collisions, Phys. Rev. D86 (2012) 054005. arXiv: 
1203.6139 doi : 10.1103/PhysRevD. 86.054005 

[11] I. Balitsky, G. A. Chirilli, Next-to-leading order evolution of 
color dipoles, Phys. Rev. D77 (2008) 014019. arXiv:0710. 
4330 doi : 10.1103/PhysRevD. 77.014019 

[12] T. Lappi, H. Mantysaari, Direct numerical solution of the co¬ 
ordinate space Balitsky-Kovchegov equation at next to leading 
order, Phys. Rev. D91 (2015) 074016. arXiv: 1502.02400 
doi:10.1103/PhysRevD.91.074016 

[13] I. Balitsky, Quark contribution to the small-x evolution of 
color dipole, Phys. Rev. D75 (2007) 014001. arXiv :hep-ph/ 
0609105 doi:10.1103/PhysRevD.75.014001 

[14] E. lancu, J. Madrigal, A. Mueller, G. Soyez, D. Triantafyl- 
lopoulos. Resumming double logarithms in the QCD evolution 
of color dipoles, Phys. Lett. B744 (2015) 293-302. arXiv: 
1502.05642, doi:10.1016/j.physletb.2015.03.068 

[15] E. lancu, J. D. Madrigal, A. H. Mueller, G. Soyez, D. N. Tri- 
antafyllopoulos, Collinearly-improved BK evolution meets the 
HERA datt arXiv: 1507.03651 

[16] I. Balitsky, G. A. Chirilli, NLO evolution of color dipoles in A = 
4 SYM, Nucl. Phys. B822 (2009) 45-87. arXiv:0903.5326 
doi:10.1016/j.nuclphysb.2009.07.003 

[17] L. Motyka, A. M. Stasto, Exact kinematics in the small x 
evolution of the color dipole and gluon cascade, Phys. Rev. 
D79 (2009) 085016. arXiv:0901.4949 doi:10.1103/ 
PhysRevD.79.085016 

[18] G. Beuf, Improving the kinematics for low-x QCD evolution 
equations in coordinate space, Phys. Rev. D89 (2014) 074039. 
arXiv: 1401.0313 doi : 10.1103/PhysRevD. 89.074039 

















